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We compare localization properties of one-dimensional Frenkel excitons with Gaussian and 
Lorentzian uncorrelated diagonal disorder. We focus on the states of the Lifshits tail, which dom- 
inate the optical response and low-temperature energy transport in molecular J-aggregates. The 
absence of exchange narrowing in chains with Lorentzian disorder is shown to manifest itself in the 
disorder scaling of the localization length distribution. Also, we show that the local exciton level 
structure of the Lifshits tail differs substantially for these two types of disorder: In addition to the 
singlets and doublets of localized states near the bare band edge, strongly resembling those found 
for Gaussian disorder, for Lorentzian disorder two other types of states are found in this energy 
region as well, namely multiplets of three or four states localized on the same chain segment and 
isolated states localized on short segments. Finally, below the Lifshits tail, Lorentzian disorder in- 
duces strongly localized exciton states, centered around low energy sites, with localization properties 
that strongly depend on energy. For Gaussian disorder with a magnitude that does not exceed the 
exciton bandwidth, the likelihood to find such very deep states is exponentially small. 

PACS numbers: 78.30.Ly 73.20.Mf 71.35.Aa; 



I. INTRODUCTION 

The term "exciton", introduced seventy- five years ago 
in the pioneering works of Frcnkeli and Wannier^ has 
become widely used to explain optical and transport 
properties of a large variety of organic and semicon- 
ductor materials^^ Within this general context, low- 
dimensional (nanoscale) systems currently attract par- 
ticular attention^ 

In low-dimensional systems, an important factor in- 
fluencing the exciton states is the presence of disorder, 
which results from static fluctuations in the host, dif- 
ferent growth conditions, as well as imperfections of the 
systems themselves. Disorder induces localization of the 
exciton states^ on certain segments of the system; the 
linear extent of such segments is usually referred to as the 
localization length. These localized states consequently 
lead to the appearance of tails in the density of states 
(DOS) outside of the bare exciton band, which are com- 
monly known as Lifshits tails^ These states play a key 
role in the low-temperature energy transport as well as 
in the optical properties of a wide spectrum of materi- 
als, such as conjugated oligomer aggregates^ and poly- 
mers^ molecular J-aggrcgates^ semiconductor quan- 
tum wells and quantum dots,— semiconductor quantum 
wires^ as well as photosynthctic light harvesting com- 
plexes^ and proteins J£ 

In discrete materials, characterized by interacting sites 
(molecular aggregates, conjugated polymers, photosyn- 
thetic complexes, spin systems), various types of disor- 
der may be considered. Commonly used are uncorrelated 
diagonal (on-site) and off-diagonal disorder, where dif- 
ferent choices still can be made for the types of disorder 
distributions; a Gaussian distribution and a box-like dis- 
tribution are the most common choices. An interesting 



alternative is provided by a Lorentzian distribution, be- 
cause the tight-binding model with diagonal uncorrelated 
Lorentzian disorder (known as the Lloyd model) is one of 
the few disorder models that allows for analytical calcu- 
lation of several physical quantities, such as the averaged 
one-particle Green's function^ and the variance of the 
Lyapunov exponent, which is a measure of the localiza- 
tion length*^ The details of the disorder model affect the 
optical response and transport properties of the above 
mentioned systems, sometimes substantially. While in 
comparison to Gaussian and box-like distributions the 
choice of Lorentzian disorder is not very common, it is 
worth noting that the latter naturally occurs in random 
systems dominated by dipolar interactions (sec Appendix 

In this paper, we perform a comparative study of un- 
correlated Gaussian and Lorentzian diagonal disorder in 
one-dimensional excitonic systems, with particular inter- 
est in the localization properties, level structure, and 
statistics of the wave functions of the exciton states in 
the Lifshits tail. We will only consider moderate disor- 
der magnitudes, where the exciton states still correspond 
to electronic excited states that are coherently shared by 
a number of molecules. There is an important differ- 
ence between Gaussian and Lorentzian disorder: the for- 
mer is characterized by a bounded second moment, while 
the second moment of the latter diverges. Distributions 
with a finite second moment give rise to exchange nar- 
rowing^ because the exciton wave functions are coher- 
ently shared by a number (N*) of monomers, they feel 
an effective disorder of magnitude V N* times smaller 
than the bare value. This effect explains the narrow- 
ness of the optical spectra of molecular J-aggregates as 
compared to their monomeric counterparts.— Exchange 
narrowing docs not occur for Lorentzian disorder^ As 
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we will show, this difference strongly affects the disor- 
der scaling of the localization properties of the excitons, 
resulting in differences in the optical and transport prop- 
erties. In addition, for Lorentzian disorder the relatively 
high density of sites with a very low energy (well outside 
of the exciton band) , also plays an important role in the 
exciton optical dynamics. 

The outline of this paper is as follows. In the next 
section, we present our model and analyze the exchange 
narrowing effect. Section UTTl contains the results of nu- 
merical simulations and discusses the local level structure 
and the statistics of the wave functions in the Lifshits tail 
for Lorentzian disorder, which we compare to previously 
obtained results for Gaussian disorder. In Sec. IIV1 we 
discuss the scaling properties of the average as well as 
the standard deviation of the localization length in the 
Lifshits tail for both types of disorder. Section [V] sum- 
marizes the paper. In Appendix lAl we show that systems 
with random dipolar interactions provide physical real- 
izations of the Lorentzian disorder model that is analyzed 
in this paper. Finally, in Appendix [Bl we present some 
mathematical details of derivations outlined in Sec. Ill Bl 



For a given disorder realization, diagonalizing the 
Hamiltonian (fl]) yields the N exciton wave functions 

\ v ) = Xm=i Vfn| n ) arL d the corresponding energies E v . 
For nonzero disorder strength, these exciton states will 
be localized on a length scale that depends on the en- 
ergy. The linear optical response is dominated by the 
states with a large transition dipole to the ground state, 

i.e. with a large oscillator strength O v = (T^—i (Pvnj , 
where we have assumed that all molecular transition 
dipolcs are equal, and the oscillator strength of a 
monomer is set to unity. These states occur in the 
neighborhood of the lower exciton band edge for the 
disorder- free system, Eb = —2 J, mostly in the Lifshits 
tail, i.e., just below E^. Their typical extension (local- 
ization length) is indicated by N* . This quantity can also 
be interpreted as the typical number of coherently bound 
molecules participating in a particular exciton state. 

There are various measures for the localization length 
of a particular exciton state. We will consider the one 
based on the inverse participation ratio 21 ' 22 i 23 or partic- 
ipation number. The latter is defined as 



II. THEORETICAL BACKGROUND 
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A. Frenkel Hamiltonian 

We consider a disordered Frenkel exciton chain of N 
molecules, described by the Hamiltonian 
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Here \n) denotes the state in which molecule n is excited, 
while all other molecules are in their ground state. E n 
denote the molecular excitation energies and — J (J > 0) 
is the nearest-neighbor interaction. We neglect interac- 
tions beyond the nearest-neighbor, as this allows for an 
analytical discussion of several important quantities. In 
Sec. lIVI we will briefly comment on extension of the model 
to include long-range dipole-dipole interactions. We ac- 
count for disorder by including a stochastic component 
in the site energies E n . Two distributions of E n will be 
considered, a Gaussian and a Lorentzian, both with zero 
mean (the transition energy of an isolated molecule is set 
to zero) 
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where and A j, denote the standard deviation and half 
width at half maximum (HWHM), respectively, which 
are measures of the disorder strength. 



which reflects the number of molecules that contribute to 
the exciton state 1^). 



B. Exchange narrowing 

A common property of systems of interacting molecules 
is that their delocalized excited states give rise to much 
narrower spectral peaks than an ensemble of noninter- 
acting molecules. This phenomenon is referred to as ex- 
change narrowing^ Its origin lies in the fact that the de- 
localized excited states do not feel the local disorder mag- 
nitude, but rather an average over its variations, which 
leads to a reduced effective disorder. More specifically, 
for an exciton state spread over a localization segment of 
length N* , the residual disorder strength is a = A/V N*, 
where A is the standard deviation (A 2 is the second mo- 
ment of the disorder distribution). 

For Lorentzian disorder, exchange narrowing does not 
occur because the second moment diverges, which results 
in the absence of the exchange narrowing effect.— Below, 
we briefly sketch these arguments. For this purpose, we 
will use the Hamiltonian JTJ) on the basis of the exciton 
wave functions of a disorder-free linear chain. 
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C. Estimates of the localization length 
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H vu i is a stochastic matrix fluctuating from one realiza- 
tion of the disorder to another. Its diagonal elements 
H vv describe fluctuations of the exciton energies due to 
disorder, while the off-diagonal part describes the scat- 
tering of excitons between different states, which eventu- 
ally results in their localization. We are interested in the 
distribution functions P{H) of these fluctuation matrix 
elements. 

For Gaussian diagonal disorder, it is given by (see Ap- 
pendix IB"]) 
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= 1/V(3/2)(7V+ 1) and k vv > = 
1/s/N + 1 {y ^ v'). Thus, the distribution 
function P vv * is also a Gaussian with stan- 
dard deviation t w A G = A G /^(3/2)(A + 1) and 
k vv i Ac = A G / VN + 1 for the diagonal and off-diagonal 
elements, respectively. The observed suppression of the 
bare disorder magnitude A G by a factor of y/N+1 
reflects the exchange narrowing effect^ 

By contrast, for Lorentzian diagonal disorder P(H) 
reads (see, again, Appendix IB"]) 
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As is seen, the distribution function of H vv i is also a 
Lorentzian. For the diagonal elements, k vv = 1, inde- 
pendently of is, i.e., P UV (H) has the same width as the 
bare distribution (|3a|) . which implies that there is no ex- 
change narrowing in this case. The off-diagonal elements 
are distributed differently, depending on v and v 1 . For 
our purpose, namely theoretically estimating the localiza- 
tion length in the neighborhood of the lower band edge, 
only k\2 is relevant (see below). The analytical result 
reads k\% = 8/(3?r), which does not show an exchange 
narrowing effect either. For arbitrary v and v' the sum 
in Eq. j8b|) can not be evaluated analytically. However, 
it can be seen that for large N it scales linearly with 
N . Thus, in the limit of large N the exchange narrowing 
effect is absent for all H,,,,i. 



The expressions in the preceding section arc valid gen- 
erally. However, only in the perturbativc limit, when 
Hvu' < \E V — E v i\, does it make sense to consider \v) 
as the (approximate) eigenstates and to interpret the 
widths of the distributions P UU {H) as linewidths for the 
absorption peak of that state. If this inequality does not 
hold, the off-diagonal matrix elements H uv i mix the ex- 
citon states, resulting in their localization on segments 
of the chain. Obviously, the perturbative limit is never 
reached for infinite chains, as then the energy separations 
of states adjacent in energy get infinitesimally small. It 
has been shown, however, that the above mixing argu- 
ments lead to an excellent estimate of the typical exci- 
ton localization length near the bare band bottom if ap- 
plied self-consistcntly to states within finite localization 
scgmentsi 24 i 25 

The self-consistent argumentation is valid for states lo- 
calized on chain segments much smaller than the chain 
length and much bigger than a single molecule (1 <C 
N* <C N) . Two notions undcrly the reasoning: (i) States 
localized on the same chain segment undergo level repul- 
sion, i.e., they have a finite energy difference. In par- 
ticular, the energy separation between the two bottom 
states in a localization segment of length N* is approx- 
imately given by the energy difference E^ — E*, where 
E* is given by Eq. ([6b)l for N = N* . Here we used the 
fact that the states resemble those of a finite homoge- 
neous chain with length N*, in particular in the sense 
that the lowest exciton state on the segment has a wave 
function without nodes, while the next higher state has 
one nodei 2 ^ (ii) Two homogeneous basis states v and v' 
localized on a chain segment of length N* are mixed by 
H* v , given by Eq. ([6c]) with N replaced by N* . The cen- 
tral argument in estimating the localization size near the 
band bottom is now that this size adjusts itself such that 
#21 = E 2 ~ E i- Namely, if < El - El, the disorder 
would only be perturbative and the states would increase 
their spread. On the other hand, if > E% — E*, the 
disorder would strongly mix the two states and would 
localize the exciton wave functions further. 

As we have seen above, for Gaussian disorder, or any 
other disorder distribution with a bounded second mo- 
ment, exchange narrowing of H21 tak es place, such that 
its typical value is ff^ = Aq/\/N*. Furthermore, as- 
suming that iV* » 1 we have E$ - El = 3n 2 J/N* 2 . 
Thus, the requirement = E^ — El leads to the esti- 
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for the typical localization size. This power-law behavior 
is in excellent agreement with previous numerical calcula- 
tionS ) 22 i 23 as well as with the analytical scaling relation 
obtained within the Coherent Potential Approximation 
(CPA)^ 



4 



For systems with Lorentzian disorder, no exchange nar- 
rowing occurs, i.e., = 8A L /(3n). Now the require- 
ment i?2i = E% — El yields 
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which reveals a different power-law scaling than for Gaus- 
sian disorder. In Sec. IIVI we will find that the A^^ 2 
scaling indeed agrees with numerical results. 

While our main interest is in the optically dominant 
band-edge states, it is interesting to apply the above ar- 
guments also to the band center, and compare to previous 
results^ Near the band center, E* ~ E*, oc 2ttJ/N* 7 
with N* now indicating the typical localization size at 
the band center. Equating this quantity to the exchange 
narrowed disorder strength, Aq/VN*, we find N* oc 
(J/ Ac) 2 . This indeed is the well known disorder scal- 
ing of the localization length in one-dimensional systems 
with Gaussian site disorder^ which in Ref. has also 
been obtained by an analytical calculation of the inverse 
participation ratio, performed within the framework of a 
one-dimensional nonlinear supermatrix a model. 



III. HIDDEN SPECTRAL STRUCTURE FOR 
LORENTZIAN DISORDER 

As we mentioned above, the optically dominant states 
in a disordered exciton chain with negative transfer inter- 
actions occur in the neighborhood of the band bottom, 
predominantly in the Lifshits tail. The exciton states 
in this tail exhibit a hidden structure,— where doublets 
of s- and p-like states often occur on the same localiza- 
tion segments (and more rarely, triplets). As a result the 
low-temperature optical response of the chain behaves 
approximately as that of a collection of segments with 
typical size N* equal to the typical localization length 
of the optically dominant states. More detailed statis- 
tics were obtained in Ref. HH, where the disorder scaling 
of the localization length, the oscillator strength, and the 
energy spacings E% — E^ were scrutinized numerically for 
Gaussian diagonal disorder. In this section, we report the 
results of similar numerical calculations with Lorentzian 
diagonal disorder and show that this both alters the hid- 
den structure near the bare band bottom and induces a 
relatively high density of strongly localized low energy 
wave functions. 

The exciton wave functions for a given realization of 
the disorder are straightforwardly obtained by diagonal- 
izing the Hamiltonian Eq. {TJ), where we pick the site 
energies E n from a Lorentzian distribution Eq. ()3a|) . In 
Fig. [TJ we show a typical realization of the wave func- 
tions in the neighborhood of the bare exciton band edge 
Eb = —2 J (upper panel) and well below it (lower panel), 
calculated for A^ = 0.05J. It is clearly seen that the 
localization properties of these two subsets of wave func- 
tions differ substantially. The states near the band edge 
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FIG. 1: A realization of the low energy wave functions near 
the band bottom (upper panel) and well below it (lower panel) 
for a chain of N = 200 sites and a Lorentzian distribution of 
site energies (HWHM Al = 0.05 J). In the upper panel three 
types of states are observed: (i) Multiplets of states (upper 
panel), (ii) states strongly localized on short segments, and 
(iii) states that form a conventional hidden level structure 
(see text for details). The solid vertical lines indicate the sites 
with large negative energy fluctuations, around which states 
deep in the tail are localized (lower panel): n = 5, n = 16, 
n = 40, n = 125, and n = 162, while the dashed vertical line 
corresponds to the high energy site n — 168. Two low energy 
exciton states (at n = 5 and n — 162) are not shown. 



are much more extended than those deep in the tail and 
exhibit a hidden multiplet structure, similar to the case of 
Gaussian disorder . 25 ' 28 The important difference is that 
for Gaussian disorder, one usually finds only singlets and 
doublets of states localized on a particular chain segment. 
In the case of a Lorentzian distribution, on the other 
hand, we see the occurrence of multiplets of three and 
even four states, often on segments with sharply defined 
boundaries. It should be stressed that this situation is 
typical; looking at various realizations, we always found 
higher order multiplets. Below, we provide an explana- 
tion for this special property of Lorentzian disorder. 

As we mentioned above, the states deep in the tail of 
the DOS (lower panel in Fig. |TJ) look different from those 
in the neighborhood of the band edge. They are, first, 
localized much more strongly (in fact, on a few sites, as 
can be seen by the eye and is confirmed by the participa- 
tion number) and they are always represented by s-like 
singlets. These states originate from a large negative 
fluctuation of one particular site energy. For moderate 
disorder magnitudes, such fluctuations occur frequently 
for Lorentzian disorder, in contrast to Gaussian disorder. 
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Such outliers in energy have another consequence: they 
act as natural barriers, providing a segmentation of the 
chain into smaller subchains. In the realization of Fig. [TJ 
these outliers occur at the positions n = 5, n = 16, 
?i = 40, n = 125, n = 162 and n = 168. Segments 
between such barriers that happen to have a length of 
around the typical localization length N* of the band- 
edge states (upper panel), can easily support the for- 
mation of multiplcts. For instance, a triplet of local- 
ized states occurs between n = 125 and n = 162, and a 
quartet is observed between n = 168 and n = 200. On 
the other hand, more strongly localized states may occur 
near the band edge as well. This happens if two closely 
spaced sites acquire large energy fluctuations, creating a 
segment considerably smaller than TV* . The states local- 
ized between n = 5 and n = 16 and between n = 162 and 
n = 168 represent two examples of this second type of 
exciton state. Finally, the segments that are appreciably 
larger than N* show a hidden level structure similar to 
that for Gaussian disorder, that is, where the localization 
segments have poorly defined boundaries and can be seen 
to overlap each other. This third type of exciton state 
can be seen in the segment between n = 40 and n =125 
in Fig. [TJ 

We proceed with a brief analysis of the wave func- 
tions of the exciton states deep in the tail of the DOS, 
which coincide with (some of) the aforementioned seg- 
ment boundaries. These strongly localized states are cen- 
tered around very low energy sites, and typically have lo- 
calization lengths, as defined by the participation ratio in 
Eq. (HJ), of around a few sites. The wave functions decay 
approximately exponentially when one moves away from 
the central low energy site. 

A simple model for these wave functions is given by a 
particle of mass m moving in a ^-function potential well. 
It is well known that such a well supports one bound state 
(with negative energy E), the wave function of which 
behaves proportionally to exp(— |x|/£) (see, e.g., Ref . I2H) . 
The relation between the penetration depth £ and the 
energy is given by2£ 



(-E) 



-1/2 



(11) 



Clearly, the meaning of £ is similar to the localization 
length N* ; both are measures of the extent of the wave 
function. The localization length calculated through the 
participation number Eq. 2] is more suited for wave func- 
tions of unknown, arbitrary shapes. The known shape 
of the wave functions for the low-energy states, however, 
suggests to consider their exponential decay lengths and 
investigate whether their energy scaling resembles Eq. 

HU 

As was shown in Refs. HO and Hi], in the neighborhood 
of the bare exciton band edge, a universal energy variable 
e exists, which is given by 
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FIG. 2: Double logarithmic plot of the numerically obtained 
penetration depth £ (in lattice units) as a function of energy 
Eq — E. The expected decay of the penetration depth with 
energy is clearly seen. The solid line shows a power law fit, 
while the dashed line represents the £ oc (Eq — i5)~ 1//2 behav- 
ior suggested by the 5-function potential model. 



where e = Ej J, £& = Eb/ J is the bare exciton band edge 
energy (ej, = —2 for nearest-neighbor interactions), and 
the disorder strength a is cither Aq / J or / J for Gaus- 
sian and Lorentzian disorder, respectively. In terms of £, 
functions of energy, like the DOS, become universal. The 
numerical coefficients a, b, and a depend on the type of 
disorder considered as well as on whether we include only 
nearest- neighbor or all dipole-dipole interactions. In par- 
ticular, for nearest-neighbor interactions and Gaussian 
diagonal disorder, this set is given by a = 0, b = 1, and 
a = 4/3(21 while in the case of nearest-neighbor inter- 
actions and Lorentzian diagonal disorder, we should use 
a = 1, b = 4 and a ~ li^L We define a reference energy 
Eq that is dependent on the disorder magnitude, which 
is the energy that corresponds to e = on the universal 
energy scale. 

An analysis of the numerically determined wave func- 
tions indeed confirms the expected general trend that 
the lower energy states are localized more strongly than 
the band edge states. For the deep tail energy states, 
the penetration depth £ was calculated by fitting an ex- 
ponential to the numerically determined wave function. 
The result is plotted in Fig. [2] against the state's en- 
ergy measured from the reference energy Eq. We ob- 
serve that these data approximately obey a power law, 
although not with quite the same exponent as suggested 
by Eq. (|TT|) . The best fit was obtained using the power 
law £(E) = 1.017 l(E — E)/J]~ 0AQ5 . The deviation of 
the numerically determined exponent from the estimated 
power law is most likely caused by the fact that the 5- 
function model assumes a continuous position variable 
and an infinitely narrow and infinitely deep well, while 
our model Hamiltonian is discrete and can thus only ap- 



6 



proximately behave as in the <5-function model. We also 
note that the very low energy states are very symmetric, 
while the exciton states that are closer to the band edge 
show more asymmetry; this makes sense, since for the 
more energetic states the surrounding sites become more 
important, at the expense of the central low-energy site. 
As the energies of the surrounding sites are random, this 
leads to the observed asymmetry. 

Summarizing this section, we point out that the na- 
ture of the localization in the Lifshits tail of the DOS for 
Lorentzian disorder clearly differs from the usual Gaus- 
sian model in a number of aspects, in particular, it al- 
ters the hidden level structure and introduces strongly 
localized states deep in the tail which are far less likely 
to occur for Gaussian stochastic variables. Close to the 
bare exciton band edge, we observe three types of states: 
strongly localized s-like states on segments much smaller 
than the localization length N*, multiplets on segments 
of lengths comparable to TV*, and finally, states that re- 
semble the conventional hidden level structure that is 
also found for Gaussian disorder. In the next section 
we provide a more detailed analysis of the scaling of the 
localization length for both types of disorder in the neigh- 
borhood of the bare exciton band edge. 



IV. LOCALIZATION LENGTH 
DISTRIBUTIONS AND SCALING 

The localization length of a given exciton state \v) can 
be calculated numerically using the participation number 
N v , Eq. It is a fluctuating quantity and is thus more 
accurately described by analyzing its probability distri- 
bution. We focus on the energy range where the optically 
dominant states reside, i.e., in the neighborhood of the 
bare exciton band edge, E\, = — 2J, and predominantly, 
just below it. To ensure a fair comparison, we fix the end- 
points of the energy interval under consideration in terms 
of the universal energy variable mentioned in Section Hm 
In the following simulations, it ranges from = —0.1 to 
e f = 0. 

The probability distribution of the localization length 
collected in this energy region, P{Ni oc ) = J2l S ( N ioc - 
N u ) (the prime restricting the summation to the selected 
energy interval), are plotted in the insets of Fig. [3] (Gaus- 
sian disorder with Aq = 0.2 J) and Fig. 2] (Lorentzian 
disorder with = 0.05 J). They show an asymmet- 
ric shape, similar to the one presented in Ref. |28| . The 
asymmetry is caused by the fact that on average, the 
states with smaller localization lengths reside deeper in 
the DOS tail, where the DOS itself is small. If we in- 
stead go towards the exciton band edge, the situation is 
reversed. 

Analyzing the moments of the various distributions 
indeed highlights the difference between the localiza- 
tion behavior for the two types of disorder, as predicted 
in Sec. HU In particular^ for Gaussian disorder, we 
find Ni oc = 5.43(A(j/J) _ ' ' , which has an exponent 
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FIG. 3: Disorder scaling of the average localization length (de- 
noted by crosses) and its standard deviation of the distribu- 
tion (denoted by the dots) for a chain of length N = 250 sites 
with Gaussian disorder. Power law fits for both have been in- 
cluded: the solid line corresponds to Ni oc = 5.43 (Ag /J) -0 ' 65 , 
and the dashed line gives SNi oc = 2.33 (Ag/ J) ~ 0,63 . Addi- 
tionally, the quotient of the standard deviation and the aver- 
age localization length is shown by the dotted line (multiplied 
by a factor of 10 for better visibility), which clearly shows this 
ratio to be constant to a good approximation. In the inset, 
we show the localization length distribution for Ag = 0.2J, 
plotted on a Ni oc scale which is normalized to Ni oc . The re- 
sulting function is universal, as it does not depend on Ag (see 
text). 



that is in good agreement with the theoretical power 
law Eq. Lorentzian disorder, on the other hand, 

leads to the scaling with a different exponent, Ni oc = 
2.19 (A^/J) -0 ' 48 , again in good agreement with the the- 
oretical estimate of Eq. (jTDJ), N* w 5.91 (A L /J) _1/2 . 
There is a deviation in the numerical prcfactors, which is 
not surprising given the arbitrariness in defining the lo- 
calization length. Most important are the different pow- 
ers found for both types of disorder, owing to the absence 
of exchange narrowing for Lorentzian disorder. 

It is remarkable that for both types of disorder, the 
standard deviation is observed to scale with almost the 
same exponent as the first moment. While the disor- 
der amplitude varies over two orders of magnitude, the 
standard deviation divided by the first moment changes 
by less than ten percent (see Figs. [3] and This im- 
plies that on the scale of the first moment the localiza- 
tion length distributions have universal shapes, which are 
presented in the insets in Figs. [3] and 2] The Lorentzian 
disorder model is observed to yield a narrower distribu- 
tion. 

All the results presented above were obtained using 
nearest- neighbor interactions. Obviously, a similar anal- 
ysis can be performed when accounting for all dipole- 
dipole interactions. Doing so induces a shift in the band 
bottom (e& = —2.404 in this case) and modifies the DOS 
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A L (J) 



FIG. 4: Same as in Fig. but now for Lorentzian dis- 
order with HWHM Al- The power law fits for the aver- 
age localization length (solid line) and the standard devia- 
tion (dashed line) shown in the figure correspond to Ni oc = 
2.19 (Al/J) -0 ' 48 and <5iV ;oc = 0.80 (Al/J) -0 46 , respectively. 
The distribution in the inset was obtained for Al ~ 0.05 J. 

in the Lifshits tail.— In turn, this changes the exponents 
in disorder scalingS ) 23 i 28 however it does not influence the 
physics that is responsible for the differences in localiza- 
tion properties of Gaussian and Lorentzian disorder. 

V. SUMMARY 

In this study, we performed a comparison of the local- 
ization properties of a one-dimensional Frenkel exciton 
model with Gaussian and Lorentzian uncorrelatcd diago- 
nal disorder, with a special focus on the energy region be- 
low the bare exciton band edge. We have found that the 
divergent second and higher moments of Lorentzian dis- 
tributions lead to a number of interesting modifications of 
the localization behavior as compared to Gaussian disor- 
der. A striking example is the absence of exchange nar- 
rowing, which, as we have shown, results in a different 
disorder scaling of the localization length from the one 
obtained for Gaussian disorder. This theoretical predic- 
tion is supported by numerical calculations, which reveal 
power law scalings with exponents that are in excellent 
agreement with the theoretical estimates. Moreover, the 
standard deviation of the localization length distribution 
scales with disorder with almost the same exponent as 
the average localization length, implying that the shape 
of the distribution is universal. 

We have also shown that the wave functions in a chain 
with Lorentzian disorder have a hidden structure that 
differs substantially from the one found for Gaussian dis- 
order. First of all, Lorentzian disorder gives rise to a rel- 
atively high density of strongly localized exciton states 
deep in the DOS tail. They resemble the bound states of 
a (5-function potential and occur as a result of large fluc- 



tuations of certain monomer transition energies; this is 
a special properety of Lorentzian distributions, resulting 
from the fact that their second and higher moments di- 
verge. These exciton wave functions are localized on only 
a few monomers; their spatial extent decreases with de- 
creasing energy The more extended states in the Lifshits 
tail, close to the band bottom, form manifolds of states 
localized on various segments of the chain. Lorentzian 
disorder produces a higher amount of multiplets as com- 
pared to Gaussian disorder, where only singlets and dou- 
blets are commonly encountered. In addition, we have 
also shown that strongly localized s-like singlet states 
may occur near the band edge on segments that are ap- 
preciably smaller than the localization length N* , while 
an exciton level structure that is comparable to the one 
for a Gaussian disorder model is found for segments 
that are considerably larger than the localization length. 
These changes may result in different behavior of the 
temperature-dependent energy transport in J-aggregates 
for these two types of disorder. 
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APPENDIX A: DIPOLAR INTERACTIONS AS A 
SOURCE OF LORENTZIAN DISORDER 

In this Appendix, we show that random dipolar inter- 
actions may result in Lorentzian disorder. Consider a 
(transition) dipole dolo surrounded by other dipoles diL 
(all of the same magnitude), where do and d\ are their 
magnitudes, while lo and L denote their orientations. We 
assume that the surrounding dipoles are randomly dis- 
tributed in a volume V and also oriented within the solid 
angle An according to a probability density /(I;), so that 
the joint probability density to find any of them some- 
where in space and somehow oriented is /(L)/V. The 
object of our interest is the probability distribution of 
the total dipole-dipole interaction of the central dipole 
dolo with the surrounding dipoles g?iL. 

M 

V^^lo.l,,!!,), (Ala) 

i=i r% 

</>(l , U, rij) = loL - 3(l n i )(l i n i ) , (Alb) 

where Vq = d^di/a 3 is a constant, is the (dimcn- 
sionless) distance between the dipoles c?q1q and diU, and 
</>(l, L, n, ) is the orientational factor with the unit vec- 
tor along Yi . The quantity V can be associated with the 
so-called solvent-induced shift of the monomer transition 
energy in molecular aggregates (see, e.g., the textbooks^ 
and @). 
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The probability distribution of V reads Furthermore, using the integral representation for the 5- 

^ function and performing the average in Eq. (|A2j) , we ob- 

P(V) = ( S ( V - £ l - n *)J J - (A 2 ) 

where the angular brackets denote the average over the 
distribution of surrounding dipoles, 

N 



n / v / dhm ~- (A3) 



j 



P(V) 



2tt 



dte vt 



<lr 
~V 



dlif(li) exp 



.Vat 



Af 



— / die 

2tt 



1-^JdrJ dkf(k) 



.Vot 

1 " ex P ( ~ l — < 



M 



(A4) 



where no = Af/V is the number density of the sur- 
rounding dipoles. In the thermodynamic limit Af — ► oo, 
V — * oo, while no = Af/V = const, Eq. (|A4[) reduces to 



1 /" c 



dte lVt 



exp 



47T TA 



(A5) 



Here, £ = dz z 2 JdnJ efl/(lj)[l — cxp(— iz0)] and 
we made the substitution z = Vot/r 3 . The integral in 
Eq. (|A5[) can be evaluated analytically. The result is a 
Lorentzian distribution 



P(V) 



1 



r 



tt (v — v ) 2 + r 2 



(A6) 



shifted from zero by Vq = (47r/3)Im£noVo and having a 
HWHM T = (47r/3)Re^n V r . Both magnitudes are de- 
termined by the dipole-dipole interaction at the average 
distance between dipoles, (47r/3)n Vo 



APPENDIX B: DERIVATIONS OF EQS. [7] AND [8] 

Here, we present the derivation of the probability dis- 
tribution for the matrix elements H vv i given by Eq. (|6c[) . 



By definition, 



P VU ,{H) = {5{H-H VV 



dte iHt { 



(Bl) 

where the angular brackets (...) = [] n=1 J dE n p{E n ).... 
denote the average over disorder realizations, with p{E n ) 
being either a Lorentzian or a Gaussian distribution func- 
tion. Further, we use for p(E n ) a representation through 
the characteristic function: 



p{E n ) 



1 

2^ 



dt 



l e iB " t « exp( 



\At n \ a ) . (B2) 



For a = 1, this formula gives a Lorentzian with HWHM 
equal to A, while for a = 2we get a Gaussian with stan- 
dard deviation A. Using this representation in Eq. (|B1[) . 
we obtain 



i r°° _ r°° i r°° 

P VV ,{H) = — dte mt T[ dt n exp(-\A a t n \ a ) — dE n e iE ^- t0 ™'^ (B3) 

27T J-oo n=l->-°° 



r 



with The integral over E n yields a (^-function, which finally 

. . 2 ttistl ixv'n , s S ivcs us 

fa* («) = Tf-T sm jTT-y sin — — . (B4) 

271" J_oo 
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where 

/ n \ V a 

fc «"" = (Zl^'WrJ ■ (B6) 

As is seen from Eq. (|B5|) . the distribution P uv i{H) is a 
Lorentzian and a Gaussian for a = 1 and a = 2, re- 
spectively, only with renormalizcd HWHM and standard 



deviation, given by k uv //S.. 

For a — 2, the summation in Eq. (|B6j) can be per- 
formed analytically for any v and i/. The result reads: 
k vv = l/ v /(3/2)(7V+l) and *w = 1/v/FTT for i/ + 
v' . In the case of a = 1, we get fc„„ = 1, while at v ^ v' 
the sum in Eq. (|F36|) depends on i/ and i/' and does not 
have a simple expression. 
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